Relaxation of quantum states under energy perturbations 



(N 



Dorje C. Brody*, Lane P. Hughston"^ and Joanna Syroka* 
*Blackett Laboratory, Imperial College, London SW7 2BZ, UK 
] Department of Mathematics, King's College London, Strand, London WC2R 2LS, UK 

(February 1, 2008) 

The energy-based stochastic extension of the Schrodinger equation is per- 
haps the simplest mathematicahy rigourous and physically plausible model for 
the reduction of the wave function. In this article we apply a new simulation 



■ methodology for the stochastic framework to analyse formulae for the dynam- 
. ics of a particle confined to a square-well potential. We consider the situation 

when the width of the well is expanded instantaneously. Through this exam- 

^ \ pie we are able to illustrate in detail how a quantum system responds to an 

' energy perturbation, and the mechanism, according to the stochastic evolu- 

■ tionary law, by which the system relaxes spontaneously into one of the stable 
. eigenstates of the Hamiltonian. We examine in particular how the expectation 

^ , value of the Hamiltonian and the probability distribution for the position of 

J> \ the particle change in time. An analytic expression for the typical timescale 

' of relaxation is derived. We also consider the small perturbation limit, and 

• discuss the relation between the stochastic framework and the quantum adi- 

59 ■ abatic theorem. 

O ■ 
(N ■ 

o ■ 

^ ; I. INTRODUCTION 

4L ' Consider a quantum system for which the Hamiltonian undergoes a sudden change h — >■ 

^ . ^ at time t = in such a manner that no mechanical 'work' is done on or by the system. 
^ ! Here h denotes the Hamiltonian before the change, and H denotes the Hamiltonian after 
the change. 

.^1 If the initial wave function of the system is given by ipo{x), the probability that the wave 

^ ■ function ^jJt{x) at a later time t > will be found in the n-th energy eigenstate Xnix) of the 

■ new Hamiltonian H is determined by the transition amplitude between ipo{x) and Xn{x)- 
The assumption that no work is done on or by the system during the sudden change in 
the Hamiltonian implies that the expectation of the Hamiltonian is preserved under such a 
transformation. That is to say, the relation 

^Q{x)h'ipo{x)dx = J ipl{x)Hil)t{x)(ix, (1) 

holds even though the Hamiltonian itself has changed. 

Just after the application of the perturbation, the wave function can be represented by 
an expansion of iPq{x) in terms of the eigenfunctions Xn{x) of H. The state then evolves 
according to the unitary law governed by H in such a way that the expectation value of H 
is preserved. According to the unitary law, the wave function iptix) will remain in a state of 
superposition of the eigenmodes of and this superposition will linger indefinitely in time. 
Indeed, it is generally necessary to append to quantum theory an additional postulate to 
the effect that only as a consequence of another kind of a sudden perturbation, namely, an 
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operation of measurement, will the wave function of the system be able to 'jump' into one 
or another of the eigenstates of the new Hamiltonian. 

In contrast, by use of the stochastic extension to the Schrodinger equation that we discuss 
in § 3, it is possible to model the dynamics of the wave function in such a manner that, after 
the system is perturbed, the wave function spontaneously relaxes to one of the eigenstates 
of the new Hamiltonian. It is a remarkable fact that the probability laws thus arising from 
the stochastic dynamics give rise to statistical predictions that are essentially in agreement 
with the standard predictions of quantum theory. 

The purpose of this paper is to present a detailed analysis of the phenomenon of stochastic 
relaxation in the simple case of a particle in a potential well when the width of the well is 
instantaneously enlarged. The structure of the paper is as follows. In § 2 we review the basic 
setup for a free particle in a potential well when the well is subjected to a sudden expansion. 
We argue that standard quantum mechanics does not give a completely satisfactory account 
of the matter. In § 3 we review the formahsm of the standard energy-based stochastic 
extension of the Schrodinger equation, and propose the use of stochastic relaxation as a 
basis for the description of the dynamics of a quantum system following a perturbation. 
In §§ 4-5 we construct the solution to the stochastic system and describe its properties, 
and in §§ 6-7 we indicate how the solution can be used for the efficient construction of 
simulations. More precisely, we show that, according to the stochastic evolutionary law, the 
system energy fluctuates randomly and eventually relaxes to one of the eigenvalues of the 
new Hamiltonian. Graphic illustrations are provided for the behaviour of the energy, as well 
as the probability density of the location of the particle in the well. In § 8 we present an 
analysis of the timescale associated with the eventuality of relaxation. This analysis then 
forms the basis of a discussion of a stochastic version of the quantum adiabatic theorem, 
presented in §§ 9-10. 



II. FREE EXPANSION IN A POTENTIAL WELL 

We analyse here the sudden expansion of a one-dimensional potential well in which a 
particle of mass /x is trapped. Before the expansion, the width of the well is L. The energy 
spectrum of the particle is given by 

where n — 1,2,... , oo, for which the corresponding eigenfunctions are 



<x<L. (3) 



We assume that initially the particle is in one of the energy eigenstates associated with the 
potential. 

Let us suppose that at t — the Hamiltonian h is changed in such a way that the width 
of the potential is increased from L to aL, where a > 1. Then for any t > the wave 
function of the system can be expressed as a superposition of the normalised stationary 
states of the new Hamiltonian H. These are given by 



2 



Xn{x) = Y ^ (S^) ' 0<x <aL, (4) 



aL \aL 
for which the associated eigenvalues are 

where n = 1, . . . , oo. It follows according to the unitary dynamics of the Schrodinger 
equation that after the expansion has taken place the system will be in an indefinite state 
of energy, and will remain so. 

If the particle is initially in the n-th eigenstate with energy e^, then the proba- 

bility Tinm that, after the expansion has taken place, the particle will be found in the m-th 
eigenstate Xm{x) of the new Hamiltonian, with energy E^, is 

'Pn{x)Xm{x)dx . (6) 

A short calculation shows that 

Clearly, if such a change of state occurs, then energy is not strictly conserved. The conser- 
vation of energy is maintained in expectation, however. Indeed, we have the identity 

oo 

^ ^ '^nm-Em (8) 
m=l 

which is valid for all n and for all a, by virtue of which we can confirm that the expectation 
value of the Hamiltonian is a constant of the motion (Bender, Brody & Meister 1999). 

To understand the energy conservation law intuitively in this example, we can regard the 
potential walls as forming the 'piston' of a one-dimensional cylinder. Suppose we consider an 
ensemble consisting of a large number of independent identical particles in the cylinder. Then 
when the piston is instantaneously moved outwards, the system undergoes a free expansion. 
During such a process, the particles do no mechanical work on the piston. As a result, 
there is no net flow of energy going out of the cylinder, and thus the energy is conserved. 
The 'energy' that is conserved in this case is the ensemble energy, i.e. the product of total 
number of particles and the expectation value of the energy of an individual particle. 

The foregoing formulae follow directly from the basic principles of quantum theory. In 
particular, the standard interpretation of quantum theory maintains that after the expansion 
the system will remain in an indefinite state of energy, until the energy is measured. However, 
there are many situations in which it is natural to presume that after the passage of some 
time the system spontaneously relaxes into one or another of the eigenstates of the new 
Hamiltonian, irrespective of whether a measurement is made. Many natural phenomena are 
of this character: after perturbation, there follows relaxation. Quantum theory, as such, 
does not account for this satisfactorily. There is arguably an implicit assumption that, in 
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the absence of specific acts of measurement, quantum systems of any significant size or 
complexity will settle into a stable eigenstate, typically an eigenstate of energy. 

In what follows we take the point of view that, after perturbation, the system does 
eventually evolve spontaneously into one of the stable eigenstates. Quantum theory in itself 
offers no clue as regards the precise mechanism according to which the system relaxes to the 
new stationary state. These questions can, nevertheless, be addressed in a very satisfactory 
way by use of the dynamics of the stochastic extension of the Schrodinger equation, as we 
shall demonstrate. 



III. ENERGY-BASED STOCHASTIC DYNAMICS 

In what follows we shall suppose that the dynamical mechanism governing the relaxation 
of the quantum system is determined by the standard energy-based stochastic extension of 
the Schrodinger equation. This is given by the following stochastic differential equation of 
the Ito type: 

d^x) = -iH^t{x)At - \a\H - Htf^t{x)dt + \a{H - Ht)^t{x)mt (9) 

for which it is assumed that there is a prescribed initial wave function il)o{x). Here Wt 
denotes a standard Wiener process, and 

^ J ip^{x)iptix)dx 

is the random process corresponding to the expectation value of the Hamiltonian operator 
H in the random state iptix)- The volatility parameter a appearing in (j^) has the units 

[a] = [Energy]-^ [Time] -^/^ (11) 

Dynamical equations of the type (||) for the evolution of the wave function and various gen- 
eralisations thereof, were introduced originally as simple models to characterise the collapse 
of the wave function when a measurement is carried out on a system (Gisin 1984, 1989, Ghi- 
rardi, et a/ 1986, 1990, Diosi 1988; see, e.g., Percival 1998, Pearle 2000 and Bassi & Ghirardi 
2002 and references cited therein for a more comprehensive account of the relevant litera- 
ture). The idea that wave function should proceed to energy eigenstates was proposed by 
Bedford and Wang (1975, 1977). The specific energy-based form of the dynamics (|^), which 
has been studied by Gisin 1989, Percival 1994, 1995, Hughston 1996, Adler & Horowitz 2000, 
Adler 2002, and Brody & Hughston 2002, amongst others, is the most parsimonious of these 
state reduction models and in many respects perhaps the most attractive as the basis for a 
fundamental model. In this paper we carry the physical application of the stochastic theory 
a step further and propose the use of (^ as an elementary model for characterising the relax- 
ation of the state of a quantum system when its Hamiltonian has been perturbed. We shall 
leave open here the question of whether the volatility parameter a governing the timescale 
of relaxation is phenomenological in character, i.e. varying according to the structure of the 
system, or universal, e.g., Planckian. Before deriving the solution to (^ and applying the 
results to perturbation theory, we shall briefly sketch some of the basic mathematical and 
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physical properties associated with the dynamics (|^). For further details, see Adler et al. 
(2001) and references cited therein. 

Let us first note that the coefficient of the Brownian motion in the third term of the right- 
hand side of @ is given by the difference of the Hamiltonian operator and its expectation, 
acting on the wave function. Thus if the system enters into an eigenstate of the Hamiltonian, 
this coefficient becomes zero, and the random fluctuations generated by the process Wt make 
no further contribution to the dynamics of ipt{x)- This property likewise applies to the second 
term, which together with the third term 'drives' the system into a state of lower energy 
uncertainty. Starting from an arbitrary initial state, the system is randomly driven into 
states with lower energy variance, until it finally reaches an eigenstate of the Hamiltonian, 
in which the variance vanishes. 

A process for which both of the coefficients of the dt term and the dWt term in the 
associated stochastic differential equation are smooth functions of the process ipti^) itself is 
called a diffusion. If the number of Brownian motions is smaller than the dimensionality 
of the process 'iptix), then the process is said to be a degenerate diffusion. In the present 
consideration, while the wave function ipt{x) is an element of an infinite dimensional Hilbert 
space, there is only a single Brownian motion driving the system, and thus the dynamical 
equation represents a highly degenerate diffusion. A typical feature of a degenerate 
diffusion, as opposed to a generic diffusion, is that it may have a 'focusing' effect. Indeed, in 
the case of (P) the wave function is focussed to one of the energy eigenstates, as a consequence 
of the dynamics. 

Given the dynamical equation @ for the wave function and the corresponding process 
(p!0|) for the expected energy, we can determine the stochastic equation satisfied by Ht. 
Specifically, this is given by 

dHt = aVtdWt, (12) 

where 

— ¥m — 

is the process associated with the variance of the energy. The variance process satisfies 

dVt = -a^^dt + a(3tdWt, (14) 

where 

mm) 

is the third central moment of the energy. As a consequence, integrating ([T^) and ([l^) we 
find that the energy process can be expressed in the form 

Ht = Ho + a [ VsdWs (16) 







and the variance process can be written as 
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(17) 



Owing to elementary properties of the stochastic integrals (|1^) and (0) one can deter- 
mine the relations 



E[if„|{if,}o<.d (18) 

and 

E[K|{K}o<sd<^t (19) 

for t < u. Here E[- ■ ■ |{Xs}o<s<t] denotes the conditional expectation given the history of 
the process Xg from time up to time t. Therefore, the conditional expectation of the 
energy process at any time u > t, given its history up to time t, is given by its value at 
time t. We thus say that Ht satisfies the martingale condition. It follows that Ht is on 
average conserved, whereas the variance Vt tends, on average, to decrease, corresponding to 
the spontaneous reduction of the system to an eigenstate. The reductive character of the 
dynamics (P) is indicated by the fact that Vt satisfies the supermartingale condition. 



IV. ANALYTIC SOLUTION TO THE STOCHASTIC EQUATION 

Despite its nonlinearity, the dynamical equation (^ can be solved exactly to yield an 
analytic solution that fully characterises the dynamical state of the system at any time in 
terms of a pair of state variables (Brody & Hughston 2002). The method of obtaining the 
solution, which is of some interest in its own right, makes use of the classical techniques of 
nonlinear filtering theory. More precisely, we consider a process of the form 

^t = atH + Bt, (20) 

where Bt is a Brownian motion, and a is the parameter that ultimately determines the 
characteristic relaxation timescale. The random variable H takes on the value Em with 
probability vr^, where vr^ is taken to be the transition probability from the initial quantum 
state to the energy eigenstate with energy Em- Intuitively, one can think of the process 
as representing the value of a phase, scaled by the constant a, together with a random noise 
term, whose strength, relative to the variable H, decreases inverse-proportionally in time. 

The random variable H is to be thought of as representing the value of the energy to 
which the system ultimately relaxes after the passage of sufficient time. Given the trajectory 
of the process C,t up to time t, we would like to determine the best estimate for the value of 
H. Because the standard error in the value of Bf grows only like the square-root of time, it 
follows that as time passes the true value of H is gradually revealed. The estimate will be 
denoted Hf, which is obtained by taking the conditional expectation of the random variable 
H, given the history of up to that time. Because C,t is a Markov process, this implies that 
Hf is the expectation of H conditional on the value C,t, that is, 

Ht = E[H\^t]. (21) 
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The argument that Ht is the best estimate for H given the history of is as follows. 
Suppose Yt is any process that at time t can be expressed as a functional of the history of 
for < s < t. Then the choice of Yt that minimises the expected mean square error 
K[{H — Yt)'^\{C,s}o<s<t] given the history of is the process defined by (0). This can be 
deduced by a straightforward variational argument. 

Because Ht is given, for each value of t, by a conditional expectation with respect to 
the random variable ^t, it follows that Ht can be expressed as a function of C,t- In order 
to determine the conditional expectation (plf ) more explicitly, we require the conditional 
probability ¥{H = Em\C,t) for the random variable H. By use of the Bayes law for conditional 
probability, we find that this is given by 

m = = V (22) 

En=l T^np{,it\H = En) 

Here p{^t\H = Em) denotes the conditional density function for the continuous random vari- 
able ^t given that H = Em- In deriving (^2]) we make use of the fact that the unconditional 
probability P(if = E^) is just t{„i for the given initial state. We have also used the relation 

oo 

p{it) = Y.''nP{m = En). (23) 

n=l 

Since Bt is a Brownian motion, it is by definition normally distributed with mean zero and 
variance t. Thus the conditional probability density for ^t is 

p{it\H = Em) = ^exp (-^(6 - crEmt)^^ . (24) 

Inserting this expression into (^), we deduce that 

, ^ TTm exp [a Emit — ^cf'^E'i.t) , , 

m = Emiit) = " rz' 1 (25) 

Since the energy process is given by the expectation 

oo 

Ht = J2EmnH = Em\it), (26) 



m=l 



it follows that 



Em=i exp {aEmit - \(y'^ElJ:) 

We note that in ( P7D the process Ht for the conditional expectation of the random variable 
H is expressed in terms of a function of t and if As a consequence, the dynamics of the 
process it can be expressed as a diffusion equation of the form 

dit = aHtdt + mt, (28) 

where Wt is a standard Brownian motion with respect to the filtration generated by the 
history of Ht- The existence of a Brownian motion Wt satisfying ( PE| ) follows as a consequence 
of well-established line of argument in nonlinear filtering theory (see Liptser & Shiryaev 
1974), the details of which in the present context are set out in Brody & Hughston (2002). 
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V. STOCHASTIC RELAXATION 



In order to see directly the relaxation of the Hamiltonian process (pTD to one of the 
energy eigenvalues, we can argue as follows. 

Let us suppose that the random variable H happens to take the value Ej. By definition 
such an event occurs with probability ttj. More precisely, we condition on the outcome of 



the event H = Ej and analyse the evolution of process (27). Then, writing = crEjt + 



we obtain, for the corresponding realisation of Ht, the expression 



_ ^jE, + Y.Z^^ ^mEm exp {a{E^ - E,)B, - \a\E^ - E,)H) 
- _ , v-oo — ——7-77^ ^TTT^ ^TT^Ix l^yj 



TC 



J + TZ^j ^-exp {a{Em - Ej)Bt - \a\E^ - Ej)H) 



where the superscript j in indicates the conditioning on the event {H = Ej}. 

Because the exponential terms appearing in both denominator and numerator have the 
property that, as t ^ 00, the probability that these terms remain positive approaches zero, 
it follows that the energy process H-l converges asymptotically to the designated eigenvalue 
Er 

Indeed for any normally distributed random variable Bt with mean zero and variance t 
it is an elementary fact that if z/ 7^ then 

lim P (exp (uBt - ^uH) > x) = (30) 

t^oo 

for any given x > 0. 

The exact, closed-form solution (^) shows how the expectation value of the system 
Hamiltonian evolves in time. In particular, it shows how the system organises itself sponta- 



neously to move into a stable, stationary state of the new Hamiltonian. Since formula (p9[) 
is expressed in terms of a simple analytic function of the Brownian motion Bt and time t, 
there arises the possibility of efficiently simulating the evolution of the process Ht. 



VI. WAVE FUNCTION DYNAMICS 

We return now to the analysis of the dynamics of the state ipt{x) of a particle trapped 
in a potential well following a sudden expansion of the width of the well. In this case the 
solution to the stochastic differential equation (RI) can be expressed in the form 



E::=i exp {-iE^t + laEUt - \cJ^Elt) x. 



Mx) = ^'^=' '^-^ — ^ ' . (31) 

(E^=i^rnexp {aEM-\o^EU))"'" 

Here Xm(a;) denotes the normalised eigenfunction of the Hamiltonian H with energy Em-, 
and the choice of the initial wave function '^q{x^ is implicit in the probability Hm- 

vTnm = ^/'o(x)xm(a;)da;^ . (32) 

The solution ( ^I]) can be verified by taking the stochastic differential of the right hand side 
of ( ^11) and using the Ito rules, and then making the substitution (^Sj) . 
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The convergence of the infinite sum in the numerator and the denominator of (^) for 
finite t may not be immediately evident. In order to check that these summations converge, 
we substitute a particular realisation of the path for ^t, say, = crEjt + Bt, when the random 
variable H happens to take the value Ej. Then on account of the fact that E^n oc m^, we find 
that the summands decay, for each fixed j, like ~ exp(— m^), and convergence is ensured. 

Let us now consider the random dynamics of the probability density function 



(33) 



for finding the particle at the location x in the interval [0,aL]. We would like to obtain 
the probability distribution for the particle in the case where the final state of the system is 
given by the eigenstate Xj{^) fo^' some given value of j. The probability for this particular 
realisation to occur is given by ttj. By rearrangement of terms, and writing 

(34) 



Em — Ej 



for the energy-level difference, we obtain 



pI{x) 



1/2 



exp (- 



mj 



t + \auJmjBt 



12 2 



t) Xr 



[X 



Em ^rn exp {(JUJmjBt 



2^ ^mj 



(35) 



for the probability density. Thus pt{x) is a measure- valued process, i.e. at each given time 
t it is a smooth density function over [0, aL], but the form of the function evolves randomly 
in time until it relaxes to the final distribution x^(x). 



VII. SIMULATION OF THE ENERGY AND THE PROBABILITY DENSITY 



A version of the paper with figures can be downloaded from 



http://theory.ic.ac.uk/~brody/DCB/sa8.pdf 



FIG. 1. Several realisations of the process Ht for the expectation value of the Hamiltonian are 
illustrated. The initial condition is chosen to be il^o{x) = (j)i{x) and the final conditions are given 
by ^poo{x) = Xi{x)^X2{x),X3{^)TX4:{^)TX5i^)TX6{x), respectively. The horizontal lines indicate the 
energy levels Ej, where we have set a = 2.5. For convenience the vertical axis is expressed in 
units of the characteristic energy e = /i^/2^L^, and the horizontal axis is expressed in units of 
the corresponding characteristic time interval When j = 5, we have vrs = 0. As a result, 

if Ej = E^ is chosen for the value of the random variable, the energy does not reduce to that 
eigenvalue. 
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With these formulae at hand, we now consider the simulation of the random trajectory for 
the expectation value of the Hamiltonian governed by the dynamics (|27|) . We also consider 
the simulation of the corresponding probability density function ( pS] ) for the position of the 
particle in the potential well. The quantum system is defined by specifying the value of 
the mass fi of the particle, the width L of the well, the volatility parameter a governing 
the stochastic dynamics, and the well expansion factor a. For simplicity, we shall examine 
the case where the system is initially in its ground state 4>i{x), although any other initial 
condition, including mixed initial states, can be treated analogously. Then the energy process 
(PPP, corresponding to the particular realisation ipooix) = Xj{^) oi the terminal state, can 
be written in the form 



Y,m T^mEm exp {(TUJmjBt - ^a^ul^-t) 

J2m '^rn exp [aUmjBt - ^a^ul^jt) 



(36) 



The associated probability density for the position of the particle in the interval [0, aL] is 
given by 



Y,m ^'^^^ exp [^auJmjBt - icr^ujl^jt) Xm{x) cos{uj^jt) 



J2m exp {aUJmjBt - \(T^Ujl^jt) 



(37) 



A version of the paper with figures can be downloaded from 



http://theory.ic.ac.uk/~brody/DCB/sa8.pdf 



FIG. 2. Several realisations of the probability distributions pt{x) for the particle in the well are 
illustrated. The initial and final conditions are V'o(a^) = (t^i{x) and t/jooix) = X2{x), respectively, 
and we set a = 2.5. The x-axis is measured in units of L, and the time axis is measured in units 
of Afi^L^/h'^a'^. Most realisations have the property that the initial state swiftly relaxes into the 
terminal state. 



The process H-l defined by (pq) can be thought of as a kind of 'Brownian bridge' that 
interpolates between the two energy levels ei and Ej. That is to say, initially the system has 
energy ei, and then it progresses along a random trajectory to reach the designated terminal 
level Ej. A similar remark applies to the density process pl{x) of (|37|), which interpolates 
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between the two functions 0f (x) and Xji^)- It should be evident that, unhke Ht, the process 
does not conserve energy. 

To obtain a better feehng for the dynamics of H-l it will be useful to examine the asso- 
ciated stochastic differential equation. In particular, if we take the stochastic differential of 
(p6|), then after a rearrangement of terms we obtain 

dHi = aV/ {E, - Hi) dt + aVMBt, (38) 

where the nonnegative process is given by 

..j _ Em T^mjEm " Hjf exp {auJmjBt - ^a'a^^.t) 
}_^f^ -Km exp [auJrnjBt - ^a^uf^jt) 

We learn from ( |38D that Hj is a mean-reverting process with mean level Ej and reversion 
rate c'^V^ . In fact, one can integrate (|38D by means of the standard technique used in the 
case of the classical Ornstein-Uhlenbeck process (see, e.g., Doob 1942) to obtain 



Hi = E, + [HI - E,) e--' y^''^-^ + a f e"'^' ""^'"'V^dB^, 

Jo 



(40) 



which shows clearly how the information of the initial condition is damped away at the rate 
a'^V^ , where 



F/ = 1 f K^d.. (41) 



t 







Alternatively, we can write (40) in the form 

Hi - E, = {hi -E,+a e'^" ^" ^''^^K'd^u) e""' ^« (42) 

which expresses the difference Hf — Ej between H^ and its terminal value as the product of 
a martingale and a positive decreasing process. 

With expressions ( p6D and (37) at hand we proceed to simulate some realisations of these 
processes. For convenience we take advantage of the fact that there is a natural energy unit 
e determined by the problem, which is given by 

(43) 



2/iL2' 

Thus, when we plot figures we can express energies in units of e, and we can express times 
in units of 

' (44) 



In the analysis that follows, energy and time will be expressed in these units. 

In Figure 1 the energy process H^ is shown for several values of j ranging from 1 to 6, 
where initially the system is in the ground state of the old potential, with energy ei = vr^. In 
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these examples the width of the potential is expanded at t = by a factor of a = 2.5. The 
most likely transition to occur, given the initial state 4>i{x), when the width is expanded by 
a factor of 2.5, is the first excited state X2{x). In this example it follows as a consequence of 
(^) that vTis = 0, and hence the transition into the fifth energy level does not occur. Instead, 
there is a transition to the fourth level. 



A version of the paper with figures can be downloaded from 



http://theory.ic.ac.uk/~brody/DCB/sa8.pdf 



FIG. 3. The expectation of i?^ is shown for a variety of realisations j = 1, 2, 3, 4, 5, 6, when the 
initial state is 4>i{x). The curious behaviour for j = 5 arises in this example because vrs = 0, and 
the state never collapses to Xbi^)- Instead, the state reduced to the eigenstate having the energy 
closest to E^, which in the present case is the fourth eigenstate. 



Several realisations of the corresponding probability distribution for the particle are 
shown in Figure 2, where the terminal wave function is the first excited state X2{x)- The 
probability of this event to be realised is as large as ~ 0.43. In most examples the initial state 
swiftly changes into the terminal state, although some interesting behaviour can occasionally 
be observed when the value of the Brownian motion Bt grows large. 

In Figure 3 the numerical average of the energy process, corresponding to the expectation 
of Hj, is shown for a range of terminal states in the case where the expansion factor is 
a = 2.5. Each graph represents the average of one thousand simulations for each energy 
level. The behaviour observed in the fifth energy level reflects the fact that, in this case, 
the terminal state is given by the fourth eigenstate X4.{^)- This is because i?4 is the closest 
energy eigenvalue to in the present example. To see this, we note that, if we write 



in the form 



12 2 



t) then as a consequence of 



it follows that can be written 



H 



^ {rCiE^ + TTmEm exp {a{um5 - ^ib)Bt - ^a' 



^ (^4 + Yl'm '^rn exp [oiuJ^^ - UJ45)Bt - ^CT^ 



m5 



UJ 



45 



li)) 



(45) 



I5 > for all 



where Yl'm ~ Sm^4- Thus, the ^-dependence cancels, and because — u, 
m 7^ 4, we see that as t — 00 the exponents in the denominator and the numerator go to 
zero and we are left with the leading term E^. 

In general, if the terminal value for the energy is chosen to be Ej when the initial state 



ipoix) is orthogonal to Xj{^)^ 



0, then the system necessarily relaxes into the 



eigenstate whose eigenvalue is the closest to Ej (S. L. Adler, private communication). 
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In Figure 4 we sketch the dynamics of the ensemble average for the probabihty density 
function for the location of the particle, in the case where the terminal state is the first 
excited state X2ix). This result corresponds to the average of a thousand examples of the 
type presented in Figure 2. 



VIII. RELAXATION TIMESCALE 

One of the advantages of the simulation methodology is that it opens up the possibility 
of a direct analysis of the timescale over which relaxation typically occurs. Let us define 
the process M^j by 



exp {auJmjBt - ^a'^uj'^jt) 



(46) 



where for simplicity of notation we suppress the time dependence of Mmj{t). The energy 
process 



conditional on the outcome Ej can then be written in the form 



(47) 



since Mjj = 1. Clearly, if M^j is sufficiently small, then Hj will approach its terminal value 
Ej and the system will have relaxed. This occurs when t > t^. 



A version of the paper with figures can be downloaded from 



http://theory.ic.ac.uk/~brody/DCB/sa8.pdf 



FIG. 4. The average of pt{x) for a thousand runs. The initial and final conditions are given by 
■0o(a;) = 01 (x) and il^ooix) = X2{x)^ respectively, and we set a = 2.5. The x-axis is measured in 
units of L, the time axis is measured in units of 4/x^L^//i^cr^, and the magnitude of pt{x) is scaled 
by a factor of 1,000. This plot illustrate how the density matrix is diagonalised. 



In order to study the scale of tr, we consider the probability that the process M^j is 
smaller than e~'^ for some number A. If uOmj > this is given by 



12 2,, 
2^ ^rnjt < 



N 



a\Umj\y/t 



mj 



a\LO, 



(48) 
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Here N{x) is the standard normal distribution function 



N{x) = ^= I e-'^y'dy. (49) 



^271 J- 



oo 



Note that in obtaining the second expression on the right-hand side of ( ^81) we have used 
the fact that Umj > 0; whereas if ujmj < we would instead have the probability 

P {Mmj < e-^) = P ( 5* > laujrnjt —] . (50) 

\ O'iOmj J 

However, owing to the symmetry relation N{—x) = 1 — N{x) satisfied by the normal distri- 
bution function, the final result (|48| ) is unaltered. 

Now, for relaxation, we would like the probability ( ^HD to be sufficiently great, say, 

P {Mmj < e-^) > 0.95, (51) 

for a suitably large A. The choice of A corresponds to how far the reduction has to proceed 
for relaxation to have effectively set in. Because N{x) ~ 0.95 when x ~ 1.65, we can put a 
bound on the time variable such that the probability (^) is greater than 0.95. This is given 

by 



la\ujmj\Vi- , ^, ^ > 1-65, (52) 

alujr, 



Tnj I 



or equivalently 



/- 3.3 + V3.32 + 8A 

In particular, for a fixed value of j, we would like this relation to hold for all values of m 7^ j. 
That is to say, we would like Mmj to become negligible for all m ^ j, which will ensure 
that — > Ej. This can be guaranteed by choosing m = j + 1, since this choice maximises 
the right hand side of (|53|) . For example, if we take A = 10, then the relaxation timescale 
obtained is of the order 

40a^ 

^ vrV2(2j + ly (^^^ 

when the terminal state is Xj{^)- Thus provided t > r^, we can be 95% confident that 
Mmj < 10~^ for every m ^ j. 

We note, incidentally, that the result in (|5^ ) is based on the assumption that the ex- 
pansion factor a is strictly and sufficiently greater than one. In the exceptional case where 
a = 1, we have tr = 0, which follows directly from (^Tj) if one notes that vr^ = for every 
m j, and thus Hf = Ej irrespective of the values of Mmj and t. 

When a ~ 1, the foregoing analysis on the relaxation timescale needs to be carried out 
more carefully. For this purpose, let us consider the case in which we have a = 1 + e, where 
the value of e > can be very small. Now, if the system is initially in the ground state 
(f)i{x), then for e ^ 1, the transition probability 7Cm{^) admits the expansion 
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4sin^(7rm) / IGsin^fvrm) 



7r^(m^ — 1)^ \7r^(m^ — 1)^ 

4(3 sin^(7rm) — 27rm sin(7rm) cos(7rm))\ . 
+ 7r2(m2 - 1)2 Je + --- (55) 

for all m. In particular, 7ri(e) is of order one, whereas the next largest transition probability 
when a ~ 1 is given, up to order e^, by vr2(e) = ^e^- The implication of this is that 
the condition for TTmMmi < to be valid for all m 7^ 1 with 95% confidence-level is 
automatically satisfied for all t > 0, and thus it is not a sufficient criteria to determine the 
relaxation timescale. 

Therefore, it is important in the small perturbation regime to determine how small e can 
be for the {p x 100)% confidence- level analysis to be viable. This can be determined if we 
replace A in (^8]) by A + \mim{e) and 0.95 in (|5T|) by p. The latter is equivalent to replacing 
1.65 in (^) by N~^{p), where N~^{x) is the inverse of the normal distribution function (^9|). 
The result gives us 



^ 2N-\p) + v/(2iV-i(p))^ + 8(A + ln7r^(e)) 



^m] I 



and we maximise the right-hand side of (^) over all m ^ j, when j = 1. The maximum is 
obtained by setting m = 2, which leads us to the timescale for relaxation into ground state. 
This is given by 

~ (n'\p) + + 41n(4e/3) + {N-^{p))') . (57) 

a UJ21 \ ' J 

In particular, for a fixed e, either A or p must be sufficient large so that the right-hand side 



of (|57|) is real. Conversely, for fixed A and p the confidence- level analysis is viable only when 



e is large enough to ensure the reality of the right-hand side of 



IX. TOWARDS THE QUANTUM ADIABATIC THEOREM 

The analysis of the relaxation timescale shows that for an infinitesimal expansion of po- 
tential well the initial eigenstate 0„(x) will almost immediately relax into the corresponding 
eigenstate Xny^\ This result suggests that when the potential well is expanded sufficiently 
slowly, then the dynamics of relaxation will force the system to remain in the n-th eigen- 
state. In other words, the stochastic evolutionary law appears to give rise to an analogue of 
quantum adiabatic theorem. In what follows we shall explore this idea in greater depth. 

For this purpose we need to consider now the dynamics of a quantum system in the 
presence of a time-dependent Hamiltonian. Let us denote by H(t) a generic time- dependent 
Hamiltonian, and for each fixed time t we write Xn{t, x) (n = 1, 2, . . . ) for the n-th eigenstate 
of the operator H{t) with eigenvalue En{t). Now, for an adiabatic approximation, we shall 
assume that H(t), Xnit,x), and En(t) are continuous in t, and vary sufficiently slowly so 
that we can write 

^^^^-iE„(t)Xn(t,x). (58) 
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Then for a general state ipt{x) the Schrodinger evolution can be written as 

= -iH{t)Mx) = -iH{t) «n(t)Xn(t, x). (59) 

n 

Therefore, we find that 

Y,('n{t)Xn{t,x) =0 (60) 

n 

in the adiabatic approximation, from which it follows at once ak{t) = 0. Thus the solution 
to the deterministic Schrodinger equation can be written as 

i^tix) = ^a„Xn(^,a;). (61) 

n 

These relations are the consequences of the adiabatic approximation ( p8|) in the deter- 
ministic unitary theory. What we would like to consider here is the implication of the 
assumption (^HD in the case of the stochastic evolutionary law (^, which in the present 
context is given by 

+\a{H{t) - Ht)i;t{^)dWu (62) 

where 

^ ^ rt{x)H{t)^t{x)dx 
* / 4)^{x)ipt{x)dx 

We note that ( |6^ ) preserves the norm of iptix)- However, for the energy process Ht we have 

dHt = Htdt + aVtdWt, (64) 

instead of (p!2|) , where 

Ht = J rt{x){dtH{t))Mx)dx. (65) 

It is natural that the process corresponding to the expectation value of the Hamiltonian is 
no longer a martingale in the time- dependent case, but rather exhibits a drift, the sign of 
which depends on the expectation value of the time derivative of the Hamiltonian. Similarly 
for the variance process we obtain 

dVt = -a'^V^dt + aptdWt 

+2 {^j rt{x){H{t) - Ht){dtH{t) - Ht)iJt{^)dx^ dt. (66) 

Now, for state reduction, a necessarily condition is that the variance process is a 
supermartingale, i.e. on average a decreasing process. We observe that the last term in (|66D 
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is given by the covariance of the two operators H{t) and dtH{t), which can be positive or 
negative. If it is negative, then the variance process Vt is a supermartingale, whereas if it 
is positive, Vt can still be a supermartingale, provided the covariance is not too large. In 
particular, if we write 



AHt = VVt (67) 

and 



AHt = \JV^T[dtH{t)] (68) 

respectively for the standard deviations of H{t) and dtH(t), then for the supermartingale 
condition we require 

2ptAHtAHt < (y\AHt)\ (69) 

where pt is the correlation between H{t) and dtH{t). If we notice the fact that — 1 < pt < 1, 
we find, to ensure that \4 is a supermartingale, it suffices that the following relation should 
hold: 

^ < lA^Htf. (70) 

In other words, we require the uncertainty in the change of H to be small compared with the 
uncertainty in H itself. Intuitively, if the change of the Hamiltonian is sufficiently slow, then 
dtH{t) will be small, and therefore we would expect AHt ^ AHt be valid in the adiabatic 
regime. In particular, relation (^) can be viewed as the consequence of adiabatic motion, 
which ensures that the variance is a supermartingale. 

Now we turn to the the consideration of the stochastic equation (^). To begin with, let 
us define the process 11^ by 

n,^ = n^exp (^aj\Ek{s) - H.,)dW, - y j\Ek{s) - Hs)Ms^ . (71) 

Then it is a straightforward exercise in Ito calculus to deduce that 

dn^ = cr{Ek{t) - Ht)U'ldWt, (72) 
from which it follows that the dynamics of the process 

«H = (nn^/' (73) 

is given by 

dakt = -\a\Ekit) - Htfattdt + \{Ek{t) - Ht)aktdWt, (74) 
Therefore, if we define 

'^t{x) =^aktXj{t^x), (75) 
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where the functions Xjit^ ^) ^^e eigenfunctions of H{t), then we obtain 



diptix) = ^ (dan) Xkit, x) + ^ ttkt (dtXkit, x)) , (76) 

k k 

and thus 

+ \a{H{t)-Ht)^t{x)dWt 

+ I 5^ an {dtXk{t, x) - H{t)xk{t, x)) J dt. (77) 



Thus far our analysis is exact. Now, in the adiabatic regime, we have the relation (^8|), 
which implies that 

{dtXk{t,x)^Hit)xk{t,x). (78) 



Therefore, substituting this relation into equation ([77|) , we see that in the adiabatic ap- 
proximation, the process ipt{x) defined by (|75|) is the solution to the stochastic dynamics 

We observe that in the adiabatic approximation the process 11^ has the interpretation 
that it represents the probability that ipt{x) is in the eigenstate Xkit, x). On the other hand, 

is also a martingale, satisfying ([72|). Thus the expected probability that ipt{x) is in the 
state Xk(t,x), given information up to time s, is precisely the probability that ipsix) is in 
the state Xk{s,x). It follows, therefore, that if ipt{x) is in the state Xfc(0,x) at time 0, it will 
be in the state Xkit,x) at time t with probability one. 



X. DISCUSSION 

Although in the present paper we have primarily examined the square-well potential, 
many of the results obtained, namely, the wave function (|3TD, the conditional Hamiltonian 
process (0), the probability distribution of the particle (|37D , and the bound on relaxation 
time (|5B|), are independent of the specific model being considered. As long as the eigen- 
values Em and the eigenfunctions Xm{x) of the Hamiltonian are known, either analytically 
or numerically, these results can be applied to study the details of the dynamics. In other 
words, the general approach outlined here, based on the nonlinear filtering methodology 
introduced in Brody & Hughston (2002), can be applied to study a wide range of problems 
in perturbation theory, including of course a more generic time-dependent Hamiltonian. We 
hope to take up this line of investigation in greater detail elsewhere. 

The point of view we have put forward in this paper is that the dynamical law (^ offers a 
simple but plausible characterisation of the subsequent evolution of a quantum system after 
the Hamiltonian has been perturbed. If a measurement of the system energy is carried out 
after a passage of time greater than the relaxation timescale, then according to the stochastic 
postulate, when the system is initially in the eigenstate 4>n{x), the eigenvalue Em will be 
observed with probability Unm- This is because the system is already in the eigenstate Xm{x) 
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with that probabihty. On the other hand, according to the unitarity postulate associated 
with the Schrodinger equation, the system is in a superposition of a myriad of eigenstates. 
Nevertheless, the outcome of the energy measurement will give with probability Unm- 
The question thus arising is whether one can distinguish the two theories by means of a 
suitable experiment. To this end, we note that, for time t > tji, the state of the system, 
according to the stochastic law, is given by a mixed-state density matrix, whose diagonal 
elements, in the energy basis, are given by vr^ = vr^m; whereas according to unitary law, 
the system is in a pure state, given by the superposition of the energy eigenstates, whose 
coefficients are given by ^/7fm. Therefore, if these pure and mixed states can be distinguished 
statistically in some way by means of experiments, then in principle one can rule out at least 
one of the two postulates indicated here. At present, the possibility of distinguishing the 
two remains an open problem. 

The authors are grateful to S. L. Adler, J. Anandan, A. Bassi, T. A. Brun, I. Buckley, and B. K. Meister 
for stimulating discussions. DCB acknowledges support from The Royal Society. LPH acknowledges the 
support and hospitality of the Institute for Advanced Study, where part of this work was carried out. 
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